Statistics of local density of states in 
the Falicov-Kimball model with local disorder 
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Statistics of the local density of states in the two-dimensional Falicov-Kimball model with local 
disorder is studied by employing the statistical dynamical mean-field theory. Within the theory 
the local density of states and its distributions are calculated through stochastic self-consistent 
equations. The most probable value of the local density of states is used to monitor the metal- 
insulator transition driven by correlation and disorder. Nonvanishing of the most probable value 
of the local density of states at the Fermi energy indicates the existence of extended states in the 
two-dimensional disordered interacting system. It is also found that the most probable value of the 
local density of states exhibits a discontinuity when the system crosses from extended states to the 
Anderson localization. A phase diagram is also presented. 
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I. INTRODUCTION 

Electron interaction and disorder strongly influence 
the properties of materials. In particular, the motion 
of charge carrier particles can be suppressed by Coulomb 
interaction and disorder, and the suppression leads to 
a metal- insulator transition (MIT). In the pure system 
without disorder the MIT can occur and is purely driven 
by electron correlations.^ The transition is commonly re- 
ferred to the Mott-Hubbard MIT. In the presence of dis- 
order, the MIT can occur even without electron interac- 
tions. The state of system changes from extended phase 
to the Anderson localization due to coherent backscat- 
tering from randomly distributed impurities.^ The Mott- 
Hubbard MIT is characterized by opening a gap of the 
density of states (DOS) at the Fermi energy, while the 
Anderson localization is a gapless insulator. At the An- 
derson localization the spectra change from continuous 
to dense discrete points. It is characterized by vanishing 
of the most probable value of the local density of states 
(LDOS),^i^ii The most probable value of the LDOS dis- 
criminates between the metal and insulator phases. The 
importance of the distribution of the LDOS in the proper 
description of the Anderson localization was stressed by 
Anderson. ^'^'^ The very distribution of the LDOS deter- 
mines the Anderson localization but not average quan- 
tities of the LDOS. From the distribution of the LDOS 
one can detect both the vanishing of the most probable 
value of LDOS as well as the gap opening of the total 
DOS, thus the distribution of the LDOS is a valuable 
tool for determining the MIT driven by both disorder 
and correlation. 

In recent years, a generalized Curie- Weiss mean-field 
theory, the dynamical mean-field theory (DMFT) was 
developed.— The DMFT essentially captures local tem- 
poral fluctuations. It has been widely applied to study 
correlated electron systems. The DMFT describes the 
Mott-Hubbard MIT wcU.'"^ However, it works with the 
arithmetic average of the LDOS, and cannot determine 



the Anderson localization in disordered systems. A sta- 
tistical variant of the DMFT, which is usually referred 
to the statistical DMFT, has been introduced to study 
systems with both disorder and interactioni^ It can 
be viewed as the DMFT formulated in real space with 
general inhomogeneous solutions.— Within the statistical 
DMFT, the self energy is a local function of frequency, 
but it also depends on the site index. In the presence of 
diagonal disorder, the self energy is also diagonal random 
variables, and gives additional dynamical random contri- 
butions to disorder. It generates a set of self-consistent 
stochastic equations. The statistical DMFT essentially 
deals with the LDOS, hence it is capable of studying 
the Anderson-Mott-Hubbard MIT. In parallel with the 
statistical DMFT, a typical medium theory (TMT) was 
also introduced to study the Anderson localization. S The 
TMT is based on the DMFT too. However, instead of 
the arithmetic average DOS, it works with the geomet- 
ric average DOS. The geometric average DOS is incor- 
porated into the self-consistent stochastic DMFT equa- 
tions, which result into self-consistent equations of the 
DMFT fashion. The TMT is essentially a mean-field 
theory of both disorder and correlation, while in the sta- 
tistical DMFT disorder is treated exactly, and only the 
correlation effects are treated in a mean-field manner. 
The TMT was employed to study the Anderson-Mott- 
Hubbard MIT in correlated electron systems with local 
disorder.— 

In the presence of disorder the LDOS forms a stochas- 
tic ensemble. The stochastic ensemble of the LDOS 
must have characteristics which discriminate between the 
metallic and insulator phases, and one has to search the 
characteristics for determining the MIT. Examples for 
the characteristics are the most probable value of the 
LDOS in the original Anderson theory of localization^!^ 
or the geometric average of the LDOS in the TMT.— How- 
ever, nearby the critical point of the MIT, quantum fluc- 
tuations may induce outliers of the statistical description 
of the stochastic ensemble of the LDOS. For a random 
sample an estimator is called robust if it is insensitive to 
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outliers." The most probable value, arithmetic average, 
median are the examples of robust estimator. The geo- 
metric average is not a robust estimator because it does 
not fulfil the linear property of the robust estimator. 
From the point of view of robust statistics the most prob- 
able value is a better estimator of a random sample than 
the geometric average. Moreover, in general, the geo- 
metric average DOS is not the most probable value of 
the LDOS, hence it does not truly represent a typical 
DOS, although it is closer to the most probable value of 
the LDOS than the arithmetic average of the LDOS. The 
geometric average DOS is sensitive to small values of the 
LDOS at individual sites, even when these values do not 
represent the most probable value of the LDOS. A decline 
of the geometric average DOS with increasing the disor- 
der strength does not necessarily imply the approach to 
the Anderson localization.— Nevertheless, when the geo- 
metric average is embedded into the self-consistent cycle 
of the DMFT, the whole method, i.e the TMT, can de- 
scribe the Anderson localizationi2ii2iiii 

The interplay between disorder and correlation in the 
MIT theory is a long standing problem. In the clean 
or noninteracting limits the MIT was well studiedi^ii^ 
However, the correlation effects in disordered systems 
still remain unclear. In particular, despite experiments 
found evidences of the metallic behavior in the two di- 
mensional electron systems, in theory it is not clear how 
electron correlations induce metallic phase in low dimen- 
sional disordered systemsJ^ In the present paper we con- 
sider the interplay of disorder and short range interac- 
tion in the MIT. Usually, the short range interaction 
is modelled by the Hubbard interaction^^ Here we take 
an alternative point of view. Rather than try to study 
the Hubbard model we take a simpler model, the spin- 
less Falicov-Kimball model (FKM).^^. The relation of the 
FKM to the Hubbard model is analogous to the relation 
between the Ising and Heisenberg models of magnetism. 
The FKM describes itinerant electrons interacting via a 
repulsive contact potential with localized electrons (or 
ions). It can also be viewed as a simplified Hubbard 
model where electrons with down spin are frozen and do 
not hop. Certainly, within the TMT the phase diagram 
of the Anderson-Mott-Hubbard MIT in the FKM and 
that in the Hubbard model share common features J^iii 
Moreover, the FKM exhibits a rich phase diagram. In 
homogeneous phase it exhibits the Mott-Hubbard type 
of MIT, although the model does not describe the Fermi 
liquid picture. At low temperature different phases with 
long-range order may exist depending on the doping and 
interaction strengthii2,ii2,i22i The FKM can also be incor- 
porated into different models to study various aspects 
of electron correlations, for instance the charge ordered 
ferromagnetism in manganites.— i^^'^'^ In the disordered 
FKM we can study different realizations of the inter- 
play of disorder and electron correlations. With local 
disorder the FKM exhibits the Anderson-Mott-Hubbard 
MITiii It will be studied in the present paper by em- 
ploying the statistical DMFT. After solving of the self- 



consistent equations of the statistical DMFT we obtain 
the distributions of the LDOS. We determine the most 
probable value of the LDOS and use it to monitor the 
Anderson localization. Wc find extended states which 
occur in the region from weak to intermediate strengths 
of interaction and disorder. For intermediate values of 
disorder or interaction there is a reentrance of the An- 
derson localization. At strong disorder the system is a 
gapless insulator, while at strong interaction the system 
is a Mott-Hubbard insulator. We also find at the cross- 
ing point from extended states to localization the most 
probable value of the LDOS exhibits a discontinuity. 

The plan of the present paper is as follows. In Sec. II 
we present the statistical DMFT through its application 
to the FKM with local disorder. Numerical results are 
presented in Sec. HI. In Sec. IV conclusion and remarks 
are presented. 

II. STATISTICAL DYNAMICAL MEAN-FIELD 
THEORY 

In this section we describe the statistical DMFT 
through its application to the FKM with local disorder. 
The Hamiltonian of the system reads 

<i,j> i i 

+Efj2fU. + uY,4cjiu (1) 

i i 

where c| (ci) and f} {fi) is the creation (annihilation) 
operator of itinerant and localized electrons at site i, re- 
spectively, tij is the hopping integral of itinerant elec- 
trons between site i and j. In the following we take into 
account only nearest neighbor hopping, i.e., tij = —t for 
nearest neighbor sites, and tij — otherwise. We will use 
t as the energy unit. U is the local interaction of itiner- 
ant and localized electrons, fj, is the chemical potential for 
itinerant electrons. It controls the electron density. Ef 
is the energy level of localized electrons. It also serves as 
the chemical potential of localized electrons and controls 
the density of localized electrons. In the following we will 
consider only the symmetric half filling case. It turns out 
that it is equivalent to — U/2 and Ef = — [//2J^ Si 
are independent random variables. They represent lo- 
cal disorder in the model. We will consider the random 
variables with uniform distribution 

^(^.)--^0(Y-k.l), (2) 

where 9(x) is the step function, and W represents the 
disorder strength. When the disorder is absent, model 
(HI) is the pure FKM. In homogeneous phase which oc- 
curs at high temperature it exhibits a Mott-Hubbard 
MYlM^^^ When the interaction is absent (f/ = 0), itin- 
erant and localized electrons are decoupled, and the itin- 
erant electron part of model ([1]) represents the Anderson 
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model, and it would exhibit the Anderson localization.-^ 
Thus when both disorder and interaction are present, 
model ([!]) would exhibit the complex Anderson-Mott- 
Hubbard MIT transition at high temperaturcr^^ The 
FKM may be realized by loading two kinds of fermion 
atoms with light and heavy masses in optical lattice. The 
light atoms play the role of itinerant electrons, while the 
heavy atoms are kept immobile as the localized electrons. 

The main idea of the statistical DMFT is to formulate 
the DMFT of the system with a realization of disorder in 
real spaccj^ When the disorder is realized, the local Green 
function is not homogeneous anymore. In this case we 
can adopt the inhomogeneous DMFT- for treating the 
interaction part. Within the approach the disorder is 
treated exactly, while the effects of electron interaction 
result into the self energy which is self-consistently calcu- 
lated by the DMFT equations. However, the self energy 
depends on the site index. The electron Green function 
can be written in real space 



(3) 



where (lo) is the self energy of the electron Green func- 
tion Gij{uj) = ((ci|c]))[^. Gq{u!) is the noninteracting 
Green function. For the FKM with a realization of dis- 



order Go(w) = [ujSij — ti 



Ei is random 



variables realized accordingly to the probability distribu- 
tion ([2]). Equation ([3]) is just the Dyson equation written 
in the matrix form. Within the statistical DMFT, the self 
energy is approximated by a local function of frequency. 
However, this local function can vary from site to site, 
i.e., 



(4) 



The approximation is strictly local. In infinite dimen- 
sions the self energy is purely local. For finite dimensions 
the approximation neglects nonlocal correlations. The 
nonlocal correlations can be systematically incorporated 
by cluster extensions of the DMFT.~ The site depen- 
dence of the self energy is generated via random vari- 
ables Ei. As a consequence the self energy Ei(w) is also 
stochastic variables. With this feature the effective mean 
field and the local Green function also are local stochas- 
tic variables. The Dyson equation ([3]) shows that the 
self energy gives additional local random contributions 
to disorder of the system. These contributions are due 
to both interaction and the interplay between interac- 
tion and disorder. However, in difference to the random 
variables Ei the contributions are dynamical. They take 
into account temporal local quantum fluctuations gen- 
erated by interaction and disorder. They also broaden 
the random energy levels generated by disorder. The self 
energy Ei(cj) is determined from an effective single site. 
Once the effective single site is solved the self energy is 
calculated by the Dyson equation 



(5) 



where Qi{uj) is the bare Green function of the effective 
single site and represents the effective mean field acting 



on site i. Gi{uj) is the electron Green function of the 
effective single site. The self-consistent condition requires 
that the Green function Gi(uj) of the effective single site 
must concise with the local Green function of the original 
lattice, i.e., 



G,{lo) = Gu{lo). 



(6) 



In the Appendix we show the exact derivation of the self 
consistent equation for inhomogeneous systems in infinite 
dimensions. Equations ©-(IHl) form the self-consistent 
system of equations for the lattice Green function and the 
self energy. They are principal equations of the statistical 
DMFT. Since the local Green function is stochastic vari- 
ables, the self-consistent equations are naturally stochas- 
tic too. In Eq. ([6]) the right hand side is a functional of 
stochastic variables G'i(w), thus the self-consistent con- 
dition ^ generates a stochastic chain of Gi{ijj) via the 
iteration process. The LDOS is defined as usually 
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Im Gri{uJ + iri), 



where = 0"*" is an infinitely small positive number. 

For the FKM the effective single site problem has the 
following action 



S, 



drdr' c\{T)g;\T - t')c,{t') 



dTU{clcJ^mT)+(3EffJf,, (7) 



where /3 = 1/T is the inverse of temperature. The parti- 
tion function corresponding to the action is 



Z,. = Tr 



(8) 



This partition function can be calculated exactly, because 
the trace over the localized electrons is independent of the 
dynamics of itinerant electrons. We obtain^S 



Zi — 2 exp 



In 



g, \iuJn) 



In 



, (9) 



where ojn — {2n + l)TrT is the Matsubara frequency. The 
Green function can directly be calculated from the par- 
tition function. Without difficulty one obtains 



Gi{iLUn) 



W, 



Oi 



+ 



li 



(10) 



where Wu = f{Ei), W^, = 1 - Wu. Here f{x) = 
l/(exp(/3a;) -f 1) is the Fermi-Dirac distribution function, 
and 



Ei — Ef 



1 



Ugi{iuJn) 



(11) 
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Note that the weight factors Woi, Wu are not simply a 
number. They are functionals of the local Green func- 
tion. One can show that Wu ~ {fjfi) is the density 
of localized electrons at site i. So far we have obtained 
the complete solution of the effective single site. This to- 
gether with the statistical DMFT equations dS])-® fully 
determine the dynamics of itinerant electrons with a fixed 
disorder realization. Within the statistical DMFT, the 
disorder is treated exactly, while the correlation effects 
are taken into account through the mean field contribu- 
tions of the DMFT. Once the self-consistent equations 
of the statistical DMFT are solved we obtain the LDOS. 
With many different realizations of disorder a data en- 
semble of the LDOS is obtained. The size of the ensemble 
depends on the lattice size and the number of disorder 
realizations. From the ensemble of the LDOS we can de- 
termine its probability distributions as well as the most 
probable value of the LDOS. The most probable value of 
the LDOS is used to monitor the MIT driven by disorder 
and correlation. 

The present approach proposes a statistical local de- 
scription of the MIT in disordered interacting systems. 
The statistical aspect is a proper description since it 
works directly with the ensemble of LDOS and use 
the most probable value of the LDOS to monitor the 
MIT.^>^>^ The local treatment in the spirit of DMFT ne- 
glects nonlocal correlations. This weakness may be seri- 
ous in low dimensional systems. However, for the two di- 
mensional FKM nonlocal correlations give nonsignificant 
contributions to the DOS in the homogeneous phasc;^^ 
The DMFT calculations for the two-dimensional FKM 
also show reasonable resultsj^ Nevertheless, the statis- 
tical DMFT can be considered as a simplest approach 
incorporating both disorder and correlation to the MIT. 



III. NUMERICAL RESULTS 

In this section we present the numerical results of the 
statistical DMFT equations. In general, the matrix in- 
version in Eq. ^ can be performed only for a finite size 
lattice. We consider a two dimensional square lattice 
with the linear size L. Thus, we perform numerical cal- 
culations for the lattice of size L and finite number 
of disorder realizations. In particular numerical calcu- 
lations were performed for L = 12 and Nd = 100. We 
take temperature T = 1 which is high enough for homo- 
geneous phase and avoiding any long-range order. Cer- 
tainly, the homogeneous phase is insensitive to temper- 
ature, however at low temperature it is unstable against 
the long-range ordered phase^^fl The Mott-Hubbard like 
MIT in the FKM occurs only for the homogeneous phase. 
Strictly speaking, the Anderson localization is well de- 
fined only at zero temperature. Here we study the in- 
terplay between the Anderson localization and the Mott- 
Hubbard MIT by investigating the single particle Green 
function at finite temperature. However, in the absence 
of long range orders the single particle Green function 




FIG. 1: The total DOS Pa{ijj) and its derivative \dpa{uj) / duj\ 
for various interactions and disorders. 



is insensitive to temperature. Hence, we can consider 
the single particle Green function at finite temperature 
as it would be at zero temperature. The symmetric half 
filling case with fixed fi = U/2 and localized electron 
density n/ = 1/2 is considered. The energy level Ef of 
localized electrons is determined accordingly by condi- 
tion n/ = J^iififi)/^'^ each, disorder reahzation. We 
solve the statistical DMFT equations in real frequency by 
iterations. The small positive number rj = 0.01 is used. 
Without disorder the FKM exhibits the Mott-Hubbard 
MIT with critical value Uc ~ 4. Thus when system is 
clean the system state is metallic for U < Uc, and is 
insulator for U > Ur- 



A. Total DOS and band edge 

First we determine the band edge of the system. Usu- 
ally, the band edge is determined from vanishing condi- 
tion of the total DOS. The total DOS is defined as the 
arithmetic average of the LDOS 



Pa(w) 



1 



(12) 



idcr 



However, strictly speaking, in the numerical calcula- 
tions the total DOS never vanishes since rj — 0.01 was 
used. In Fig. [T] we plot the total DOS and its deriva- 
tive \dpa{u!) /du!\. One can see that at the band edge the 
total DOS sharply changes and its derivative exhibits a 
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FIG. 2: (Color online) Probability distribution P{pi) of the 
LDOS Pi at the Fermi energy for various disorder W in the 
weak interaction case {U = 2). 



pronounced peak. We use the position of the peak to 
determine the band edge. The value of the total DOS 
at the band edge, phe, also serves as the cutoff of the 
DOS. Below this value p^e the DOS approximately van- 
ishes. For strong interaction and weak disorder, the total 
DOS opens a gap at the Fermi energy. It is similar to 
the Mott-Hubbard insulator in the clean case. For such 
cases we also use the peaks of \dpa{uj) / duj\ to determine 
the band gap. For strong interaction and strong disor- 
der, the gap opened at the Fermi energy closes. As we 
will see later, the system is still localized, but gapless. It 
is a crossover from the Mott-Hubbard to the Anderson 
insulator by disorder. 



B. Anderson- Mott-Hubbard MIT 

The probability distribution of the LDOS is con- 
structed from statistical data of the LDOS which is ob- 
tained after solving the statistical DMFT equations. In 
Fig. [2] we present the histogram of the probability dis- 
tribution of the LDOS at the Fermi energy for [/ — 2 
and various disorders. This value of interaction (U = 2) 
corresponds to the metallic phase in the clean limit. 
For weak disorders the probability distribution has a 
monomodal structure with a shaped peak at its most 
probable value. In this regime the most probable value 
of the LDOS has a finite value, thus the system is in 
an extended state. This is an unambiguous evidence of 
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FIG. 3: The most probable value pmpv, arithmetic average pa 
and geometric average pg of the LDOS at the Fermi energy 
via disorder strength W in the weak interaction case ((7 = 2). 



the existence of extended states in the two-dimensional 
disordered interacting system. As the disorder strength 
increases the peak broadens, and then a second peak is 
developed. Thus the probability distribution has a bi- 
modal structure. The high value mode is nearly fixed 
almost independently on the disorder strength. The low 
value mode moves towards to zero value. Actually, as the 
disorder strength increases the low value mode becomes 
the most probable value of the LDOS, and it approxi- 
mately vanishes at some value of disorder. The vanishing 
of the most probable value of the LDOS is detected in the 
sense that its value is below the density cutoff pbe- The 
vanishing of the most probable value of the LDOS mani- 
fests the Anderson localization. Thus the system exhibits 
a MIT from extended state to the Anderson localization 
as the disorder strength increases. Perhaps, the bimodal 
structure of the probability distribution of the LDOS is 
due to special features of the FKM. The low value mode is 
due to the Anderson localization when disorder increases. 
The high value mode reflects the non- Fermi liquid behav- 
ior of the FKM in the weak interaction regime. For weak 
interactions the chemical potential remains pinned at the 
effective level of / localized electrons. In the presence of 
disorder most of the LDOS still persist with the pinning 
property, thus most of the LDOS at the chemical poten- 
tial keep the same value. The bimodal structure may be 
absent in the systems where the Fermi liquid properties 
are maintained. 

In principle, we can determine the most probable value 
as the value at which the probability distribution reaches 
its global maximum. However, the probability distribu- 
tion is constructed by histogram, its most probable value 
is sensitive to the width of histogram bars. To avoid the 
inaccuracy, we determine the most probable value by the 
half sample mode algorithmi^^ This algorithm is a fast 
routine for locating the most probable value of a finite 
statistical sample. The half sample mode algorithm is 
based on finding the smallest interval that contains half 
number of the sample points. The most probable value 
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FIG. 4: (Color online) Probability distribution P{pi) of the 
LDOS Pi at the Fermi energy for various disorder W in the 
intermediate interaction case {U = 6). 



must lie in the obtained half sample. Repeat this half 
sample procedure until obtain the half sample with two 
or three sample points. Then one can easily locate the 
most probable value of the sample. In Fig. [3] we plot the 
most probable value as well as the arithmetic and geo- 
metric average of the LDOS at the Fermi energy for com- 
parison. It shows as the disorder strength increases the 
most probable value shifts from the higher value mode to 
lower value one. After a critical value of disorder strength 
(Wc « 11.1) the most probable value approximately van- 
ishes, thus the system changes to the Anderson localized 
phase. At the crossing point the most probable value ex- 
hibits a discontinuity. The most probable value of the 
LDOS is not an order parameter of MIT, since it does 
not associate with any symmetry breaking, but the dis- 
continuity of the most probable value of the LDOS may 
be considered as a sign of the first order phase transi- 
tion. Figure [3] also shows that both the arithmetic and 
geometric average of the LDOS never coincide with the 
most probable value of the LDOS. The arithmetic and ge- 
ometric averages of the LDOS monotonously decrease as 
the disorder strength increases. However, the decreases 
do not necessarily imply the approach to the Anderson 
localizationii^ 

In Fig. U we present the histogram for the probability 
distribution of the LDOS at the Fermi energy for various 
disorders and U — 6. This value of interaction {U = 6) 
corresponds to the insulator phase in the clean limit. In 
difference to the weak interaction case, at weak disorders 
the probability distribution of the LDOS at the Fermi en- 



ergy has almost a delta function like structure. It means 
that almost all LDOS at the Fermi energy vanish. The 
total DOS opens a gap at the Fermi energy. It is anal- 
ogous to the Mott-Hubbard insulator in the clean limit. 
As the disorder strength increases, the delta peak broad- 
ens, and the probability distribution of the LDOS shows 
a long tail. However, the most probable value of the 
LDOS still approximately vanishes, while the arithmetic 
and geometric average are finite. The total DOS now 
closes the opened gap at the Fermi energy. The system 
state is still localized, however gapless. It is analogous to 
the Anderson localized phase in the noninteracting limit. 
We still refer it to the Anderson localization, although 
the physics nature may be different. In general, at finite 
interaction and finite disorder there are no precise def- 
initions of the Mott-Hubbard and Anderson insulating 
phases.— These phases rigorously exist only in the clean 
or noninteracting limits. In disordered interacting sys- 
tems, both phases are characterized by vanishing of the 
most probable value of the LDOS at the Fermi energy. 
However, the total DOS of the Mott-Hubbard insulator 
opens a gap at the Fermi energy, while the one of the 
Anderson localization is gapless. The scenario of clos- 
ing the opened gap by disorder in the strong interaction 
case can be understood from the atomic limit. In the 
clean atomic limit each site has two energy levels, one 
single occupancy and one double occupancy, separated 
by interaction strength U. At the half filling, each site 
is occupied either by one itinerant electron or by one 
localized electron. The double occupancy levels remain 
empty, thus the charge gap is equal to U. When disorder 
is added, each of these two energy levels is shifted by ran- 
domly fluctuating site energy —W/2 < Si < W/2. For 
W < U the situation remains unchanged, thus there is 
still a charge gap for electron excitations. When W > U, 
the double occupancy levels at some sites may be shifted 
lower than the single occupancy levels. Thus the double 
occupancy levels of a fraction of sites are either occu- 
pied or empty. As a result the charge gap is closed. The 
phase may be interpreted as a mixture of Anderson and 
Mott-Hubbard insulators. 

In contrast to the weak interaction case, where the 
probability distribution of the LDOS has the bimodal 
structure, in the intermediate and strong interaction 
cases the probability distribution of the LDOS keeps its 
monomodal structure. As the disorder strength increases 
further, the most probable value of the LDOS first in- 
creases from zero value and then decreases back to zero 
value. In the first stage the system changes from the lo- 
calized state to extended state, while in the second stage 
the system changes back from the extended state to lo- 
calized state. It is a reentry effect of the Anderson local- 
ization. In Fig. [5] we plot the most probable value, the 
arithmetic and geometric average of the LDOS at the 
Fermi energy as a function of the disorder strength for 
[/ = 6. In difference to the weak interaction case, as the 
disorder strength increases, the most probable value and 
the averages of the LDOS increase from zero value, reach 
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FIG. 5: The most probable value Pmpv, arithmetic average pa 
and geometric average pg of the LDOS at the Fermi energy via 
disorder strength W in the strong interaction case (i7 = 6). 
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FIG. 6: Phase diagram for the Mott-Hubbard- Anderson MIT. 
The dotted line separates the phase region with finite gap in 
the total DOS. 



their maximum, and then decrease. It also shows that the 
geometric average of the LDOS approximately vanishes 
only in the Mott-Hubbard insulator phase. However in 
this phase the arithmetic average and the most proba- 
ble value of the LDOS approximately vanish too. In the 
Anderson localized phase only the most probable value 
of the LDOS vanishes. In the region of intermediate val- 
ues of the disorder strength, the most probable value of 
the LDOS is finite. This evidence unambiguously indi- 
cates the existence of extended states for intermediate 
disorders. It also shows that disorder can drive the sys- 
tem from the insulating state to extended one. However, 
the insulating state should be the gapless localized state. 
Disorder cannot drive a Mott-Hubbard insulator directly 
to extended state. One may speculate the MIT scenario 
of the intermediate interaction case as a screening of dis- 
order which leads to close the gap at the Fermi energy, 
and then the standard scenario of the Anderson localiza- 
tion as in the weak interaction case. At the transition 
point the most probable value also shows a discontinuity 
like in the weak interaction case. We can use the dis- 
continuity to detect the crossing point from extended to 
localized states. 

In Fig. [6] we plot the phase diagram. It clearly dis- 
tinguishes three phase regions. Extended phase is those 
states that the most probable value of the LDOS at the 
Fermi energy is finite. The insulator phase is charac- 
terized by the vanishing of the the most probable value 
of the LDOS at the Fermi energy. This phase is sep- 
arated into the Anderson localization where the total 
DOS is gapless and the Mott-Hubbard insulator where 
the total DOS opens a gap at the Fermi energy. The 
Anderson localization occurs for strong disorder, while 
the Mott-Hubbard insulator occurs for strong correla- 
tion. The extended phase appears only for weak and 
intermediate disorder and correlation. For a weak dis- 
order, as the interaction increases, the system changes 
from the extended phase to the Anderson localization, 
and finally to the Mott-Hubbard insulator. For an in- 



termediate value of disorder, as the interaction increases, 
the system changes from the Anderson localization to the 
extended phase, and then back again to the Anderson lo- 
calization. This is a reentry effect of the Anderson local- 
ization. For a weak interaction as the disorder strength 
increases the system changes from extended to localized 
phase. For intermediate and strong interactions there is 
a crossover from the Mott-Hubbard insulator to the An- 
derson insulator by closing the gap at the Fermi energy 
by disorder. For a fixed intermediate interaction there is 
also the reentry effect of the Anderson localization when 
the disorder strength is varied. The phase diagram is 
qualitatively analogous to the one calculated within the 
TMT.i^ Although the geometric mean alone does not in- 
dicate the Anderson localization, its fully embedding in 
the self consistent cycle of the DMFT may describe the 
Anderson localizatioui ^i^^i^^ However, within the statis- 
tical DMFT there is a discontinuity of the most proba- 
ble value of the LDOS at the phase boundary, whereas 
within the TMT the geometric average is continuous at 
the phase boundary. The two limiting cases W — and 
U = are special. In the clean limit W — 0, there is 
the Mott-Hubbard MIT, although the metallic phase is 
not a Fermi liquid. The noninteracting limit U — is 
controversial in two dimension. Our result agrees well 
with the real space renormalization group calculations f2E 
However, scaling theory did not find any true metallic be- 
havior in the system. '^'^ There is only a crossover from ex- 
ponentially to logarithmically localized states. Certainly, 
the present phase diagram is constructed from statistics 
of the LDOS, and the actual transport properties still 
remain unclear. Nevertheless, it was demonstrated that 
the two dimensional disordered Hubbard model can have 
a delocalizing effect ^i^^i^i^i 
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FIG. 7: The statistical average of the most probable value 
(mpv) and the arithmetic average (a) of the LDOS at the 
Fermi energy via the number of disorder realizations. The 
error bars are their standard deviations. {U — 2, L — 12, 
iVi, = 5). 
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FIG. 8: The most probable value pmpv and the arithmetic 
average pa of the LDOS at the Fermi energy via the lattice 
size L. The number of disorder realizations for each lattice 
size is given in Table I. {U = 2). 



C. Finite size effects 

The results of the previous subsections basically per- 
mit finite size effects. There are two sources of the finite 
size effects. One is the finite size L of the lattice, and 
the second is the finite number Nd of disorder realiza- 
tions. First, we study the finite size effects of Nd- We 
fix the lattice size L — 12, and consider different num- 
bers of disorder realizations. For each Nd we generate 
Nb different bins of size V^Nd for disorder realizations. 
For each bin we calculate the most probable value as well 
as the arithmetic average of the LDOS at the Fermi en- 
ergy. Then we calculate their statistical average p„ and 
standard deviation ctq, i.e., 

n— 1 

n— 1 

where a denotes the most probable value (mpv) or the 

arithmetic average (a), pa"' is the most probable value 
or the arithmetic average of the LDOS, obtained from 
the nth bin calculations. In Fig. [7] we plot the statistical 
average and the standard deviation of the most proba- 



L 


8 


10 


12 


14 


16 




240 


150 


110 


80 


60 


N 


15680 
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TABLE L The lattice size L, the number of disorder realiza- 
tions Nd, and the total number A*', which are used for study 
of the finite size effects. 



ble value or the arithmetic average of the LDOS at the 
Fermi energy for various Nd and Nj, = 5. We choose 
two values of disorder nearby the transition point. One 
is in the extended phase, and the other is in Anderson 
localized phase. Figure [7] shows the arithmetic average 
of the LDOS (the total DOS) has little fluctuations even 
for Nd — 50. It almost independent on Nd already from 
Nd — 50. The most probable value of the LDOS fluctu- 
ates in the extended phase more strongly than in the lo- 
calized phase. Certainly, in the localized phase the most 
probable value of the LDOS approximately vanishes, that 
the fluctuations of the vanishing value is negligible when 
Nd varies. The most probable value of the LDOS seems 
to be reasonable from Nd = 100. Thus at Nd = 100 the 
finite size effects of Nd are small and do not significantly 
change the results. 

Next, we study the finite size effects of the lattice size. 
The DMFT calculations for clean systems show the fi- 
nite site effects are small and controllable.® In disordered 
systems, one can notice that the finite size effects of the 
ensemble of the LDOS depend mostly on the size of the 
ensemble, i.e., on the number N = L^Nd- N is the total 
number of LDOS which are obtained in numerical cal- 
culations. Therefore to study the finite size effects of 
L alone, when the lattice size L is varied, we have to 
change Nd accordingly, that the total number N keeps 
more or less the same value. In Table I we present sev- 
eral values of L, and corresponding values of Nd that the 
total number N is around 15000. We use the parame- 
ter values in Table I for study of the finite size effects of 
L. For all cases the whole bin of disorder realizations is 
kept more or less the same. We also choose two values 
of disorder nearby the transition point. One is in the 
extended phase, and the other is in the localized phase. 
The most probable value and the arithmetic average of 
the LDOS at the Fermi energy via the lattice size are 
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plotted in Fig. [H It shows that the arithmetic average 
of the LDOS is ahiiost independent on the lattice size, 
at least from L — 8. The most probable value of the 
LDOS also slightly fluctuates as the lattice size varies. 
These fluctuations mainly are due to the statistical fluc- 
tuations of flnite size of disorder realizations. As we al- 
ready showed in Fig. [7l the statistical fluctuations of the 
most probable value of the LDOS in the extended phase 
is larger than in the localized phase. This feature is con- 
sistent with Fig. [3 where the most probable value of the 
LDOS fluctuates in the extended phase stronger than in 
the localized phase. Both results show that the flnite size 
effects in our study are small and do not change signifi- 
cantly the picture of MIT. It is interesting to note that 
the DMFT can be performed for flnite size lattices, and 
the obtained results are not significantly different from 
the ones of the thermodynamical limitj^ 



IV. CONCLUSIONS 

In this paper we have studied the Mott-Hubbard- 
Anderson MIT in the two-dimensional FKM with local 
disorder by the statistical DMFT. Within the statisti- 
cal DMFT the correlation effects are resulted into addi- 
tional dynamical local random variables, which are self- 
consistently determined from the local single site dynam- 
ics. The probability distribution and the most probable 
value of the LDOS are calculated. The localized phase 
is detected by vanishing condition of the most probable 
value of the LDOS at the Fermi energy. The scenario 
of the MIT in the system depends on the interaction. 
For weak interactions, which correspond to the metallic 
phase in the clean limit, the system changes from ex- 
tended to localized states as the disorder strength in- 
creases. For intermediate interactions, which correspond 
to the insulating phase in the clean limit, as the disor- 
der strength increases the system crosses from the Mott- 
Hubbard to the Anderson insulator, and then it transits 
to extended states and goes back again to the Anderson 
localized phase. Thus there is a reentrance of the An- 
derson localization. At the crossing point from extended 
to localized states the most probable value of the LDOS 
exhibits a discontinuity. For strong interactions only lo- 
calized states exist. There is only a crossover from the 
Mott-Hubbard to the Anderson insulator by closing the 
opened gap at the Fermi energy by disorder. The results 
also confirm the delocalizing effect in the two dimensional 
disordered interacting system. However, the phase dia- 
gram was determined only by statistics of the LDOS, and 
the actual transport properties of the system still remain 
unclear. We leave the problem for further study. 
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APPENDIX A: SELF CONSISTENT EQUATION 
OF THE INHOMOGENEOUS DMFT IN 
INFINITE DIMENSIONS 

In this Appendix we present the derivation of the self 
consistent equation of the inhomogeneous DMFT in infi- 
nite dimensions. Formally, we can derive the self consis- 
tent equation for inhomogeneous systems in the same way 
as for homogeneous systems^^ Certainly, the derivation 
for homogeneous systems is based on the cavity method 
and the Hilbert transformj^. Since the cavity method for 
homogeneous systems is also formulated in real space, we 
can follow it closely. First a disorder realization is fixed, 
and then all fcrmions are traced out except for a single 
site I. In the infinite dimensions we obtain the Green 
function which represents the effective mean field^ 

G^^iii^n) = iLUn + fJ- - ei ~ Ai{iu!n), (Al) 

Ai{iuj„) = J2tuG'l^^{tUn)tji, (A2) 

where g[^^ (itUn) is the Green function of the model with 
site / removed. A;(ia;„) can be considered as a hybridiza- 
tion function. The cavity Green function can be ex- 
pressed through the original lattice Green function^ 

Gy(iw„) = Gy(iw„) \ . (A3) 

Inserting Eq. (|A3p into Eq. (|A2p one obtains 

Ai{iuJn) = [t ■ G{iuJn) ■ t]j^ 

[t.G(»c.„)]„[G(»o;„)-t]„ ^ (^4) 

Gu{iuJn) 

where t is the hopping matrix. The lattice Green function 
can be rewritten as 

G{iLJn) = [l(it^n) - t - ^{iUJn)]~\ (A5) 

where {icun) — {iujn + M — £i)Sij . In infinite dimensions 
the self energy is purely local)^ hence the self energy ma- 
trix S(ia;„) is diagonal. For homogeneous systems the 
terms in the right hand side of Eq. (|A4p are calculated 
by using the Fourier and Hilbert transforms.^ For inho- 
mogeneous systems the Fourier and Hilbert transforms 
are replaced by the matrix multiplication and inversion 
in real space. One can notice that 

t-G-t t-G-[^-S - G-^] 
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= t • G • - S) - t, (A6) 
t • G = - S - G-^] ■ G 

= (^-S).G-1. (A7) 

Using relations (|A6p - (|A7p . from Eq. (|A4p we obtain 

+ (l(«^n) - 5](zcj„)) • G{iujn) ■ (|(icj„) - S(iw„)) 
- [(|(iw„) - S(iw„)) • G(iw„) - 

[G{iuj^) ■ {^{lujn) - " l]JGn{iuJn). (A8) 

Here we have used tu — 0. Since both matrices ^(iwn) 
and 'S{ioJn) are diagonal, from Eqs. (|Aip and (jASP we 
obtain 



Equation (|A9p is just the self consistent equations (O-®. 
It is exact in infinite dimensions for any inhomogeneous 
system. The above derivation of the self consistent equa- 
tion can be considered as a generalization of the homo- 
geneous one. In infinite dimensions, the sum in Eq. (|A2p 
runs over infinitely many neighbouring sites, that the hy- 
bridization function is replaced by its average value. As 
a consequence, all spatial fluctuations of the environment 
surrounding a given site are suppressed, thus the Ander- 
son localization is prohibited.— At infinite dimensions the 
Anderson localization is absent in any disordered inter- 
acting system. 



Gl ^{iuJn) = T,iXiuJn) + l/Gu{iuJn)- 



(A9) 
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